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SIMULATION OF JET QUENCHING IN HEAVY ION COLLISIONS 
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The method to simulate the rescattering and energy loss of hard partons in heavy ion collisions has 
been developed. The model is a fast Monte-Carlo tool introduced to modify a standard PYTHIA 
jet event. The full heavy ion event is obtained as a superposition of a soft hydro-type state and 
hard multi-jets. The model is applied to analysis of the jet quenching pattern at RHIC and LHC. 
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1. Introduction 

One of the important tools for studying the 
properties of quark-gluon plasma (QGP) in 
ultrarelativistic heavy ion collisions is the 
analysis of a QCD jet production. The 
medium-induced energy loss of energetic par- 
tons, "jet quenching" , should be very dif- 
ferent in the cold nuclear matter and QGP, 
resulting in many observable phenomena . 
Recent RHIC data on high-p-r particle pro- 
duction at y/s = 200 A GeV are in agree- 
ment with the jet quenching hypothesis (see, 
e.g., 2 and references therein). At LHC, 
a new regime of heavy ion physics will be 
reached at ^/snn = 5.5^4 TeV where hard 
and semi-hard particle production can stand 
out against the underlying soft events. The 
initial gluon densities in Pb+Pb reactions at 
LHC are expected to be much higher than 
those at RHIC, implying a stronger partonic 
energy loss, observable in new channels. 

In the most of available Monte-Carlo 
heavy ion event generators the medium- 
induced partonic rescattering and energy loss 
are either ignored or implemented insuffi- 
ciently. Thus, in order to analyze RHIC data 
on high-pT hadron production and test the 
sensitivity of LHC observables to the QGP 
formation, the development of adequate and 
fast Monte-Carlo tool to simulate the jet 
quenching is necessary. 



2. Physics model and simulation 
procedure 

The detailed description of physics model can 
be found in our recent paper 3 . The approach 
bases on an accumulating energy loss, the 
gluon radiation being associated with each 
parton scattering in the expanding medium 
and includes the interference effect using the 
modified radiation spectrum dE/dl as a func- 
tion of decreasing temperature T. The basic 
kinetic integral equation for the energy loss 
AE as a function of initial energy E and path 
length L has the form 
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where I is the current transverse coordinate 
of a parton, dP/dl is the scattering proba- 
bility density, dE/dl is the energy loss per 
unit length, A is in-medium mean free path. 
The collisional loss in high-momentum trans- 
fer limit and radiative loss in BDMS ap- 
proximation 4 (with "dead-cone" generaliza- 
tion of the radiation spectrum for heavy 
quarks 5 ) are using. We consider realistic 
nuclear geometry 6 and treat the medium in 
nuclear overlapping zone as a boost-invariant 
longitudinally expanding quark-gluon fluid. 
The model parameters are the initial con- 
ditions for the QGP formation for central 
Au+Au (Pb+Pb) collisions at RHIC (LHC): 
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the proper formation time r and the tem- 
perature To. For non-central collisions we 
suggest the proportionality of the initial en- 
ergy density £o to the ratio of nuclear over- 
lap function and transverse area of nuclear 
overlapping 6 . The simple Gaussian parame- 
terization of gluon angular distribution over 
the emission angle 9 with the typical angle of 
the coherent radiation 0q ~ 5° 7 . 

The model has been constructed as 
the Monte-Carlo event generator PYQUEN 
(PYthia QUENched) and is available via In- 
ternet 8 . The routine is implemented as a 
modification of the standard PYTHIA_6.4 
jet event 9 . The event-by-event simulation 
procedure includes the generation of the ini- 
tial parton spectra with PYTHIA and pro- 
duction vertexes at given impact parameter, 
rescattering-by-rescattering simulation of the 
parton path length in a dense zone, radia- 
tive and collisional energy loss per rescatter- 
ing, final hadronization with the Lund string 
model for hard partons and in-medium emit- 
ted gluons. 

The full heavy ion event is simulated as 
a superposition of soft hydro-type state and 
hard multi-jets. The simple approximation 3 
of hadronic liquid at "freeze-out" stage has 
been used to treat soft part of the event. 
Then the hard part of the event includes 
PYQUEN multi-jets generated according to 
the binomial distribution. The mean num- 
ber of jets produced in AA events at a given 
impact parameter is a product of the num- 
ber of binary NN sub-collisions and the inte- 
gral cross section of hard process in pp col- 
lisions with the minimal transverse momen- 
tum transfer p™ m . The extended in such a 
way model has been also constructed as the 
fast Monte-Carlo event generator 10 . Note 
that ideologically similar approximation has 
been developed in . 



3. Jet quenching at RHIC 

In order to demonstrate the efficiency of 
the model, the jet quenching pattern in 
Au+Au collisions at RHIC was considered. 
The PHOBOS data on ^-spectra of charged 
hadrons 12 have been analyzed at first to 
fix the particle density in the mid-rapidity 
and the maximum longitudinal flow rapidity, 
F max = 3 5 (fi gure l). To calculate multi- 
jet production cross section, the higher or- 
der corrections factor K = 2 was used. The 
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Fig. 1. The pseudorapidity distribution of charged 
hadrons in Au+Au collisions for three centrality sets. 
The points are PHOBOS data, histograms are the 
model calculations. 
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Fig. 2. The transverse momentum distribution of 
neutral pions in Au+Au collisions for three centrality 
sets. The points are PHENIX data, histograms are 
the model calculations. 
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rest of the model parameters have been ob- 
tained by fitting PHENIX data on pr-spectra 
of neutral pions 13 (figure 2): the kinetic 
freeze-out temperature Tf = 100 MeV, max- 
imum transverse flow rapidity Y^ nax = 1.25 
and p^ in = 2.8 GeV/c. The nuclear modifi- 
cation of the hardest domain of px-spectrum 
was used to extract initial QGP conditions: 
T = 500 MeV and r = 0.4 fm/c. Figure 3 
shows that our model reproduce pt~ and cen- 
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Fig. 3. The nuclear modification factor Raa for 
neutral pions in Au+Au collisions for two central- 
ity sets. The points are PHENIX data, histograms 
are the model calculations. 
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Fig. 4. The azimuthal two-particle correlation func- 
tion for pp and for central Au+Au collisions. The 
points are STAR data, dashed and solid histograms 
are the model calculations for pp and Au+Au events 
respectively. 



trality dependences of nuclear modification 
factor Raa (determined as a ratio of parti- 
cle yields in AA and pp collisions normalized 
on the number of binary NN sub-collisions) 
quite well. 

Another important tool to verify jet 
quenching is two-particle azimuthal correla- 
tion function C'(Aip) - the distribution over 
an azimuthal angle of high-py hadrons in the 
event with 2 GeV/c < px < p^' s relative to 
that for the hardest "trigger" particle with 
p£ ig > 4 GeV/c. Figure 4 presents C(Atp) 
in pp and in central Au+Au collisions (data 
from STAR 14 ). Clear peaks in pp collisions 
at Aip = and Aip = ir indicate a typi- 
cal dijet event topology. However, for cen- 
tral Au+Au collisions the peak near tt disap- 
pears, ft can be interpreted as the observa- 
tion of monojet events due to the absorption 
of one of the jets in a dense medium. Fig- 
ure 4 demonstrates that measured suppres- 
sion of azimuthal back-to-back correlations 
is well reproduced by our model. 

We leave beyond the scope of this paper 
the analysis of such important RHIC observ- 
ables as the azimuthal anisotropy and par- 
ticle ratios, which are sensitive to the soft 
physics. In order to study them, a more care- 
ful treatment of low-py particle production 
than our simple approach is needed (the de- 
tailed description of space-time structure of 
freeze-out region, resonance decays, etc.). 

4. Jet quenching at LHC 

The developed model was applied to analyze 
various novel features of jet quenching at the 
LHC (see 15 < 16 < 17 for details). Let us just to 
enumerate some main issues. 

Jets tagged by leading hadrons. The 
relation between in-medium softening jet 
fragmentation function (measured in lead- 
ing hadron channel) and suppression of jet 
rates due to energy loss out of jet cone, was 
analyzed in our work 15 . We have found 
that the specific anti-correlation between two 
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effects allows one to probe parton energy 
loss mechanism (small- angular radiation vs. 
wide-angular radiation and collisional loss). 

Jets induced by heavy quarks. The possi- 
bility to observe the medium-modified frag- 
mentation of hard b-quarks tagged by a en- 
ergetic muon in heavy ion collisions was an- 
alyzed in the work 16 . We have found that 
the reasonable statistics, ~ 10 4 events per 
one month of LHC run with lead beams, can 
be expected for the realistic geometrical ac- 
ceptance and kinematic cuts. The significant 
softening b-jet fragmentation function deter- 
mined by the absolute value of partonic en- 
ergy loss and the angular radiation spectrum 
is predicted. 

■Z/y* +jet production. The channel with 
dimuon tagged jet production in heavy ion 
collisions was analyzed in the work 17 . The 
correlations between pair and jet, as 

well as between [i + fi~ pair and a leading 
particle in a jet, were first numerically stud- 
ied. We have found that the medium-induced 
partonic energy loss can result in significant 
Pr-imbalance between pair and a lead- 

ing particle in a jet, which is quite visible 
even for moderate loss. 

5. Conclusions 

The method to simulate jet quenching in 
heavy ion collisions has been developed. The 
model is the fast Monte-Carlo tool imple- 
mented to modify a standard PYTHIA jet 
event. The full heavy ion event is obtained 
as a superposition of a soft hydro-type state 
and hard multi-jets. The model is capable of 
reproducing main features of the jet quench- 
ing pattern at RHIC (the pr dependence 
of the nuclear modification factor and the 
suppression of azimuthal back-to-back cor- 
relations). The model was also applied to 
probe jet quenching in various new channels 
at LHC energy: jets tagged by leading par- 
ticles, b-jets, dilepton-jet correlations. The 
further development of the model focusing on 



a more detailed description of low-pr particle 
production is planned for the future. 
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